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Abstract 

We introduce a new class of nonlinear equations admitting a repre- 
sentation in terms of Darboux-covariant compatibility conditions. Their 
special cases are, in particular, (i) the "general" von Neumann equation 
ip — [H,f(p)}, with [f(p),p] = 0, (ii) its generalization involving certain 
functions f(p) which are non-Abelian in the sense that [f{p),p] 7^ 0, and 
(iii) the Nahm equations. 

1 Introduction 

An investigation of collective phenomena in quantum mechanics leads to various 
nonlinear evolution equations. Nonlinear equation of a Schrodinger type was 
derived as a phenomenological equation for the order parameter in superfluid 
HC4 [0, ||. Recent experiments on Bose-Einstein condensation || significantly 
raise an interest in nonlinear generalizations of the Schrodinger equation (for a 
review see [Q). Another kind of nonlinear Schrodinger equations, a by-product 
of work on classification of groups of diffeomorphisms || , was recently related 
to certain aspects of D-brane dynamics || . 

In more realistic situations, where entanglement between interacting parti- 
cles is properly taken into account, one does not arrive at nonlinear Schrodinger 
wave equations but rather at their density matrix (von Neumann-type) nonlin- 
ear versions [Q 

-iX=[X,h(X)]. (1) 
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In this case the Hamiltonian h(X) is considered as a "non-Abelian function" of 
the density operator X. As is well known, equations analogous to (|l|) are often 
encountered in quantum optics and field theory if one deals with the Heisenberg- 
picture evolution of observables. 

Still another class of nonlinear von Neumann type equations may be de- 
rived in dissipative contexts M, fif or on the basis of various entropic variational 



principles llG, 11 



Of some interest is the fact that for a special class of h(X) Eq. (|l|) can be 
rewritten as 

iX=[H,f(X)}. (2) 

Nonlinear equations of this type appeared in the frameworks of nonlinear 
Nambu-type theories (l2) and nonextensive statistics |l3) . It should be stressed 
that f(X) does not always take the usual form known from spectral theory |Q. 
We will refer to the equations of the general form (Q) , (|J) as nonlinear equations 
of the von Neumann type. Eq. (|J) acquires an additional fundamental flavor if 
one recalls that for X = \ip)(ip\ and for all functions constructed via the spec- 
tral theorem, which satisfy /(0) = and /(l) = 1, one finds f(X) = X and 
therefore the dynamics of X is equivalent to the linear Schrodinger equation. 

Similar nonlinear equations can be found also in classical theories. The best 
known physical example is the Euler equation for a freely rotating rigid body 

X = [H,X 2 ]. (3) 

The more abstract versions are related to the Euler-Arnold equation for an 



"iV-dimensional rigid body" 15 , the Lie-Poisson equations occurring in fluid 
dynamics 15, lq], and the iV-wave equations for electromagnetic waves in non- 
linear media ||L7| . 

Particularly interesting and in recent years very intensively investigated class 
of nonlinear equations are the Nahm equations Jl8| . Their solutions are used as 
an intermediate step in construction of non-Abelian monopoles. One family of 
solutions is in a one-to-one relationship to the Euler rigid-body equations. In 
this sense the Nahm equations may be regarded as a kind of generalized von 
Neumann equations. 

Finally, quite recently nonlinear equations on free associative algebras, in- 
cluding the ones of the form (|l|), (|J) with h(X) and f(X) being (noncommuta- 
tive) polynomials, were considered in the framework of the symmetry approach 
to classification of integrable ordinary differential equations Jigj] . It was found 
in particular that Eq. (||), where f(X) = iX 3 , is symmetry for Eq. (||). A class 
of equations discussed in was termed "non-C-integrable" . Below we show 
that some of them belong to our class of integrable equations with Darboux 
covariant Lax representation. It is worth mentioning that we do not assume a 
polynomial structure of the RHS of Eqs. ([I]), (||). 
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|2]| , |22| is an appropriate 



27 1 or rather of its analogue 



As we can see, the reasons for generalizations of the linear von Neumann 
equation may be different, but they all finally lead to the same fundamental 
difficulty: The resulting equations involve a large (often infinite) number of 
degrees of freedom and effective integration procedures are difficult to find. The 
situation is additionally complicated by constraints typical of density matrices 
or Hamiltonians. 

It was only recently that soliton methods were applied to the density ma- 
trix version of (||) J|(], ^l], |2^] . The progress was made possible by the obser- 
vation that there exist Darboux-covariant Lax representations of certain von 
Neumann-type equations. The technique used in 
modification of the dressing method |23|, |||, |26| , 
constructed via a binary Darboux transformation 



lim, Hi- The technique is 
called the Darboux transformation since the construction of generalized gauge 
transformation is performed with the help of additional solutions of the Lax 
pairs. 

The Darboux-type method of integrating the density-matrix analogue of 
Eq. (||), introduced in |^(| and further generalized in pi| , led to discovery of 
the so-called self-scattering solutions. The process of self-scattering continuously 
interpolates between two asymptotically linear evolutions. This very character- 
istic property was found in all the nontrivial solutions of nonlinear von Neumann 
equation obtained by the above technique. 

The paper presents further development of our results. All the equations we 
derive can be regarded as compatibility conditions for Darboux-covariant Lax 
pairs. Previously discussed nonlinear von Neumann equations as well as the 
Nahm equations are particular examples of the class under consideration, but 
they form just a tip of an iceberg. 

Darboux covariance of Lax pairs is proved in detail along the lines of pl| . An 
alternative proof, involving a more general class of Darboux transformations, 
is given in this volume in |B2]. The two constructions are based on different 
mathematical techniques, show different aspects of the same problem, and it is 
not completely clear whether they are entirely equivalent. 

The layout of the paper is as follows. The compatibility-condition repre- 
sentation of a family of nonlinear equations of the von Neumann type is given 
in Sec. II. The compatibility conditions are brought into a closed form by a 
special choice of operator coefficients of the Lax pair. The coefficients are de- 
fined in terms of additional functions satisfying restrictions following from the 
compatibility conditions. A wide class of such functions is proposed. Examples 
of nonlinear equations of the von Neumann type that are generated by some of 
these functions are given in Sec. III. Darboux covariance of the Lax pairs with 
the operator-valued coefficients is proved in Sec. IV. In the next section we show 
that the restrictions carrying the compatibility condition into integrablc nonlin- 
ear von Neumann equations are Darboux-covariant if the functions introduced 
in Sec. II are used. 
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2 Darboux-integrable equations 



We begin with the overdetermined system of linear equations (the Lax pair) 

-ij> = ^a{\) . . 

where A and z\ are complex numbers, ip is an element of a linear space L, 
A(X) and -ff(A) are linear operators L L belonging to an associative ring, 
the dot denotes a derivative (i.e. an operator satisfying the Leibnitz rule). The 
compatibility condition for the Lax pair is 

iH(X) = [A(X),H(X)]. (5) 

Assume the operators entering the Lax pair are rational functions of A with 
operator coefficients 

l M 

A(X) = + ( 6 ) 

fc=0 k=l 
N 

H(X) = "£X k H k . (7) 

k=0 

The compatibility condition implies two sets of relations between operators B k , 
C' k and H k 

N 

[H k ,B m - k ] = (N <m<L + N), (8) 

fc-max{0,m-L} 
min{Ar,m+Af} 

53 [ff fc ,C fc _ m ] - (-M<m<0) (9) 

k=0 

and the system of differential equations 

rn min{N,m+M} 

— iH m = 53 [Hk>Bm-k]+ 5Z [Hk>Ck-m] (10) 

fc-max{0, m — L} k—m+1 

for (0 < m < N). In order to reduce Eqs. ( |l0| ) to equations of the von Neumann 
type one needs to write them in a closed form. In general, Eqs. (||), (||) are 
inconvenient for defining operators B k and C k in terms of H k . Nevertheless, 
one can express B k and C k explicitly through operator H k by imposing on them 
some additional relations which obey Eqs. (||), (|J). It is clear that not all such 
additional relations have to be consistent with the requirement of Darboux- 
covariance of the Lax pair. 
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Consider 



B k = 
Ck 



1 



-jL — k 



{L-k)\ V* 

(. 



~L-k 



c=0 



jM-k 



(M-jfe)! \de M - k 



(11) 
(12) 



where f(X,X) and g(X, A) are properly denned functions of operator X and 
parameter A. The operator at the RHS of the first equation of the Lax pair now 
reads 



fc=0 v 

M ~ 1 \k-M / Ak \ 

fe=0 v 



c=0 



(13) 



There exists a large class of functions f(X, A) and g(X, A) that results in opera- 
tors Bk and Ck identically satisfying conditions (^) and ([)]). The class is defined 

by 



[/(X(A),A),X(A)] = [ 5 (X(A),A),X(A)] = 0. 



(14) 



To prove the covariance of Eqs. ( p.2[ ) under the binary Darboux trans- 

formation we also assume that these functions possess an additional property, 
namely they are covariant with respect to the similarity transformation: 



f{TXT-\ A) = Tf(X, X)T-\ g(TXT-\ A) = Tg(X, X)T~\ 



(15) 



where T is a transformation. The above conditions are satisfied, for example, 
by polynomials in X and sums of negative powers of polynomials in X . If X is 
selfadjoint, the same is valid for all f(X) defined via the spectral theorem. 

In such a case Eqs. (0) and ([)]) turn out to be identically fulfilled as a 
consequence of the trivial identities 



jjL[f(H(X),X),H(X)] 



A=0 



0, ^[g(H(X),X),H(X)] 







and Eq. dlG) can be written in equivalent form 



N 



iH m = J2 [Hk,B„^ k }+J2i H ^ C k- m } (0<m<N). 



(16) 



k—m+l 



fc=0 



In the next section we will see that there are two representations of the com- 
patibility condition and they correspond to Eqs. (Q) and (||). 
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3 Examples 



Below we present a few examples of integrable nonlinear von Neumann-type 
equations that correspond to different choices of positive integers N, L, M and 
functions f(X, A), g(X, A). In what follows we will use the notation 



If N = 1 Eqs. (On) imply 



Hx =H, H a = p. 



H = 0. 



3.1 

N = 1, L = 1, f(X, A) = I"(n6 N), A) = 
The compatibility condition gives the equation 

'n-l 



ip 



Y j H n - k - 1 pH k ,p 



k=0 



(17) 



The Darboux-covariant Lax pair for this equation was found in pl| ] . For n = 2 
Eq. (Il7|) reads 



ip=[pH + Hp,p}= [H lP 2 ] 



(18) 



which is equivalent to Eq. (j^). Mutual replacement of -ff(A) = p + \H and 
A{\) = Hp+pH+XH 2 in the corresponding Lax pair results in the compatibility 
condition 

i(Hp + pH) = [H,p% 



which is essentially a form of Euler's top equations given in JTa, MS 



3.2 

N - 1, f(X, A) = 0, M = 1, g(X, A) = g(X). 
Here we have 

*P = [9(p),H] . 

The Lax-pair representation and Darboux covariance properties of this equation 
have already been established in It should be stressed that the function 

g(X) is basically arbitrary. The cases g{p) — ip 3 and g{p) — ip^ 1 were consid- 
ered in H|. 



G 



3.3 

N = 1, L = 3, f(X, A) = \- 2 (a X 2 + (b + Xbx)X 3 + (c + Aci + A 2 c 2 )A 4 ), 
g(X,A) = 0. 

The compatibility condition becomes 

ip = [h(p),p] = [H,F(p)}, (19) 

where 

h(p) = a {pH + Hp)+b {pH 2 + HpH + H 2 p)+b l {p 2 H + pHp + Hp 2 ) 
+c (pH 3 + HpH 2 + H 2 pH + H 3 p) 
+ Cl (p 2 H 2 + pHpH + pH 2 p + Hp 2 H + HpHp + H 2 p 2 ) 
+c 2 (p 3 H + p 2 Hp + pHp 2 + Hp 3 ), 
F(p) = a p 2 + b (p 2 H + pHp + Hp 2 )+b 1 p 3 

+c Q (p 2 H 2 + pHpH + P H 2 p + Hp 2 H + HpHp + H 2 p 2 ) 
+d(p 3 H + p 2 Hp + pHp 2 +Hp 3 )+ c 2 /A 

Here ao, bo, b\, cq, c\, c 2 are arbitrary complex parameters independent of A. If 
the dot is a derivative with respect to a time variable t, they can depend on t. 
The same is also valid for the next example. Let us note that the map p i— > h(p) 
is not a function of p in the standard sense of the spectral theory Jl4j (such as 
g{p) of the previous subsection). In particular, [h(p),p] ^ 0. We refer to such 
maps as non-Abelian functions, or non-Abelian nonlinearities . 

3.4 

N = l, f(X,X) = 0, M = 2, 

g(X, A) = (a + Aai) ((bo + A6i)l + xy 1 + (c + Ac x ) ((d + Adi)l + X) _1 ((e + 
Xejl+X)' 1 . 

In this case we obtain 

ip = [H,F(p)}, (20) 

where 

F{p) = aoibvl + py^hl + H^bol + p)- 1 -a^bvl + p)- 1 
+c ((d l + p) _1 (dil + H)(d l + p)-\e Q l + p)' 1 
+ (d l + py\e Q l + p)- x (eil + H)(e l + p)' 1 ) 
-ci(d l + p)" 1 (eol + /3)" 1 . (21) 

As opposed to the previous examples F(p) is a non-Abelian nonpolynomial 
function. 
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For N = 1 the nonlinear equations involve only two types of operators: p 
and H. Increasing N we can introduce non-Abelian nonlinearities involving an 
arbitrary number of different operators. 

3.5 

N = 2, L = 2, f(X, A) = X 2 , g(X, A) = 0. 

This is the simplest example of TV = 2 nonlinearity. The compatibility 
conditions are 

ip = [H 2 ,p] + [H 2 ,p 2 } = [H 2 + H 2 p + pH 2 ,p], 
iH = [H 2 ,Hp + pH], 

H 2 = 0. (22) 

This system is equivalent to a nonlinear von Neumann equation with two types 
of nonlinearity: One given in an implicit form and the other of the Eulcr type. 



3.6 

N = 2,L = 1, f(X, A) = X, g(X, A) = 0. 
The Lax pair is 



zxil>x = ^(^o + AFi+A 2 ^), (23) 
-ij>\ = Va(#i + AF 2 ) (24) 



with the compatibility conditions 



H 2 = 0, (25) 
m = [H it H ], (26) 
iH Q = [H lt H ]. (27) 



Defining 



F\ - (H -H 2 )/(2i), (28) 

F 2 = {H + H 2 )/2, (29) 

^3 = H 1 /(2i) (30) 

and the connection V/ = / + [/, F3] we can write the compatibility conditions 

as 

VFi = i[F 2 ,F 3 ], (31) 

VF 2 = i^a.-Fi], (32) 

VF 3 - i[F u F 2 ]. (33) 
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The connection can be trivialized if we find an invertible solution £ of the linear 
problem 

£ = -£fs. (34) 

Then 

fk = WC 1 (35) 
satisfies the standard Nahm equations 

A = i[/ 2 ,/ 3 ], (36) 

A = i[/ 3 ,/i], (37) 

/a = (38) 



4 Binary Darboux transformation 

The first step towards extending the technique of Darboux transformations to 
integrable nonlinear von Neumann-type equations on associative rings is to es- 
tablish the Darboux covariance of the Lax pair without any additional con- 
straints. In this section we show that an appropriate formulation of the binary 
Darboux transformation makes the Lax pair covariant. 

Assume x is a solution of the Lax pair with parameter v. 

| -ix = xA{v) 

and ip is a solution of the dual Lax pair with parameter /x: 

( = ■ (40) 

We further suppose that these systems can be related with an operator P sat- 
isfying P 2 — P and 

-iP = PA(v)P ± - P±A{fi)P, (41) 

where P± = 1 — P. The above assumptions are fulfilled, for example, if X = (x\ 
and ip = \ip) are some "bra" and "ket" associated with a Hilbert space. Then 

P = MM 
<x\v>) ' 

Another example is provided bymxn matrix X a- n d n x m matrix ip. P is then 
defined by 

p = i p{x ( p)' 1 x- 
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Some realizations of Darboux transformations in infinite dimensional cases were 
given in 34 . Very recently a new construction of the Darboux transformation 
in terms of Clifford numbers was described in p5| . 
Denning 

ff,[l] = ipD\, (42) 

v — u , 

D x = 1 + y P 43 

fi — A 

we come to the following 

Theorem 1. The Lax pair (^) with the coefficients defined by Eqs. (Q), (^) is 
covariant with respect to the binary Darboux transformation {ip, A{\) 1 H(\)} — > 
{^[l],A(A)[l],if(A)[l]}, where 

L M 

A(X)[1] = ^A fc i3 fe [l]+^^C fc [l], (44) 

fc=0 k=l 
L 

B k [l] = B k + (n-p) ]T (v m - k - 1 P±B m P-v m - k - 1 PB m Px)(*5) 

M 

C k [l] = C k - (n -v)Y (n k - m - l P±C m P - v k - m - x PC m P^) (46) 



m—k 



and 



N 



H{X)[1] = (47) 



k=Q 



H k [l] = H k 



N 



■(ji-v) J2 (l* m - k - 1 PxH m P-v m - k - 1 PH m P x ) (48) 



m—k+l 

Proof: The condition of covariance of the second equation of the Lax pair with 
respect to the transformation yields 

H(X)[1] = D- 1 H(X)D X (49) 

N N N 

' v — A u — A 

k=0 k=0 r» k=0 

Taking into account 

PH(v)P ± = P±H(jjl)P = 
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we are able to rewrite the previous expression in the following manner 

N N 

ff(A)[l] = E^ + ^E^ -" k )PH k P± 

fc=0 k=0 
N 

+qv(A t -, 1 fc )w 

li — A — 

^ k=0 

N N k-1 

fc=0 k=l 3=0 

N k-1 

+0* - v ) E E A* _J '~ V'PxfffcP, (50) 

fc=l 3=0 

which is equivalent to Eq. ([47|). 

From the condition of the Darboux covariance of the first equation of the 
Lax pair we have 

A(X)[1] ^ D- 1 A(X)D X - iD-'Dx. 
Substitution of Eqs. §l|), (||) gives 

Mm = ( 1+ ^p)a(a)(i + ^p 

(l + (PA(^)Px - P_ AMP) 

= A(A) + ^P(A(X) - A(^))Pj_ + ^Z£p x (A(X) A(ji))P 

v — X [A — X 

L M 

= E^ + E>* 

fc=0 k=l 

(L k-1 M k-1 \ 

2 E ^ J '"V'Pfc -EE ^-^-"Ck p 
fe=i i=o fc=i j=o y 

(L fe-1 A/ fe-1 \ 

fc=l 3=0 fc=l j=0 / 

The final expression is (|44|)-B 

5 The main theorem 

Theorem 1 establishes Darboux covariance of Lax pairs involving operators with 
positive or negative powers of spectral parameters. The compatibility condition 
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is also covariant: Transformed operators Cfc[l] and Hk[l] solve Eqs. (||), 

([}[), and ([To|). However, if there are additional relations between the operators, 
they do not have to be Darboux covariant. 



Theorem 2. The relations ( |ll| ) and ( |12[ ) are Darboux covariant if Eqs. ( p.5| 
are fulfilled. 

Proof: Let us check the Darboux covariance of Eq. (pdj), i.e. 



3k[l] 



1 



■:/. In'. V.'/s' 
It follows immediately that 



5=0 



Bl[1] =B L 



and 



5=0 



5=0 



k\{v — [L)V 



fc-l ; 



For k ^ L we have, using Eq. ( }49| ) 
1 



(£-fe)! V* L_fe 



1 



(L-jfe)! V* L_fc 



jL-fc 



5=0 



/ d 



^ L — k a 



(L-fe)! 



jL — k—c 



-d: 



L — /c a 1 

EE 



5=0 v* 



d° 



-b 



—f^H^^) 



a=Q f>_-0 
L — k a 

EE 

a=l 6=0 



(L-k-a)\b\ \ds L - k - a 



jL — k—a 



-D7} 



B 



d b 

5=0 V* 
0* 



6^A 



L+6-a 



jL — k—c 



-D: 



Br 



5=0 



i •"" " ^ d^ Dl/ " 



5=0 



5 = 



5 = 



+ (f — A*)^ 



L-fc-l 



PS, 



L — fe a 

= EE 



6-1 



a=l 6=1 v ' v 



jL—k—a 



-d: 



Bh+b-aP 



5=0 



12 



L-fc 



V 1 / d L-h-a 

+ 2s (L-k-a)\ l* L - fc - a V< 

a— 1 v 



B L - a + {u- /i)i/ L - fc - 1 PB L 



L — fe— 1 a 

EE 



(/i - *V 



6—1 / r ]L — k—a 



^ (L - k - a)\ \dc;L-k-a D lh 
a=l 6=1 v y x 

L-fc 

+ ^( A1 -^)/- 1 i? fe+b P 

| / d L-k-a 



BL+b-aP 



C=0 



6=1 
L-fe-1 

- — Z) -1 

^ (L-fc- a)! U<r L - fe - a V« 

a=l x ' v 



5 



L— a 



?=0 



+ S fe + (i/ - n)v J - k - 1 PB L 
= E X> - *V"> - n)v L - k -^PB L+b _ a P 

a=l 6=1 
L-fe 



+ ^( M -^ 6 - 1 Sfe+6P 
6=1 

L-fe-1 

+ ]T {v- V i) V L - k - a - 1 PB L ^ a ^B k ^{v- V i)v L - k - 1 PB L 



<L-k 



L-k-1 



Sfe + (M-^) E^ lflfe + bP - E V L - k - a - 1 PB L -a 



,6=1 



a=0 



L-fe-1 



a=l 6=1 / 
'L-k L-k-1 \ 

Sfe + (M - f) ( E ^ 1 P±B k+b P - v L - k - a - 1 PB L -aP±+Sk) , 



where 



,6=1 



L-fe 



a=0 



L-fe-1 



L-k-a-lp BL a p 



4 = ^/-^Sfe+tP- X! " 

6=1 a=0 

L — £;— 1 a 



a=l 6=1 



L-fe 



L-fe-1 



= j2» b ~ lpBk +>> p - E " 

6=1 a=0 

L — £;— 1 a 

+ J2^ 1,yL ^ aPB L+b-aP 



L-k-a-lp B p 



a=l 6=1 
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L—k—1 a 



a=l 6=1 



Combining, respectively, the first and the third, the second and the fourth terms 
gives 



L-k 



L-k-l 



a=l b=l 

Finally, we obtain 



a=0 b=0 



1 



(L — k)\ V* i_fc 



-ih — k 



'L-k 



?=0 
L-k-l 



,L — k — rn—l 



PB L . m P x (51) 



\m= 1 



m=0 



The last expression coincides with Eq. (pis]). 

Let us prove the Darboux covariance of Eq. (|lj) : 



1 



(. 



M-k 



ft[11 =(M^M ff(£)[ll£) 
One can show that 



e=0 



d k 

de k 
d k 



D, 



de k 



d: 



6=0 



6=0 



5 kQ l + k\(y - v)ti- k - l P, 



(52) 



Then 



1 / d M ~ k 



(M-k)l \de M - k - 



6=0 



1 



( d 



M-k 



g{D- 1 H(e)D £ ,e) 

(D- 1 g(H(e),e)D, 
(M-k)l 



6=0 



(M-k)\ \de M - k 
1 / d M - k 
(M-k)\ \de M - k 

_^ M-k a 

(M - k)\ ^ ^ \M -k- a)\al (a - &)!&! 



6=0 



( d 



M—k—a 



\ de M-k-a 



d: 



6=0 



de b 



D, 



6=0 
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M—k a 



= ^ ^ (M - k - a)M 



a=0 6=0 



x (V-fc-a)ol + (M — k — a)!( M - v)v k+a - M - la p) 
x C M+ b-a{s bo l + b\{v - M)M" 6-1 ^) 



M-fc 

[l — V 



+ - /.) E l^ 1 1 + — P Cn. fc P 

6=1 \ v J 

M-fe-1 , 

a=0 V M / 

M-fe-1 a 



a=l 6=1 



= (1 + ^P W 1 + ^P 



v J \ /i 

(M-fe-1 M-fe \ 

E v k+a ~ M - l PCM- a P±_ - E /i- 6 - 1 i , ±Cn. fc P + A fc , 
a=0 6=1 / 

where 

M-fe-1 M-fe 

A fc = ^ ^ +0 - M M - 1 PCM-aP - E ^"^PCb+kP 
o=0 6=1 
M-fe-1 a 

+ (""/*) E E^^'^^^^^-aP 
a=l 6=1 

M-fe-1 M-fe 

= E ^"^V^M-aP - E t^^PCb+kP 

a=0 6=1 
M-fc-1 a 

+ E E^" b_li/fe+a " Mpc ^-« p 

a=l 6=1 

M-fe-1 a 



o=l 6=1 



Combining, respectively, the first and the third, the second and the fourth terms 
we obtain 



M-fe-1 



A * = E E^ b ~ v+a ~ Mpc M + 6- a p 



a=0 6=0 
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M—k a 



J2 J2^ b » k+a ~ M ~ lpC M +b -aP = 0. 



a=l 6=1 

Finally, we have 

1 / d M - k 



(M-k)\ \de M - k 



9(H(e)[l],e) 



/M-k-l M-k \ 

(PL -u)[ v k+a - M - l PC M - a P± - £ ^Pj.C^P 



6=1 



C k + ^PC k P ± + v -^p 1 _c k p 




= C k + (p - v) [J2 iy k+a - M - 1 PC M ^ a P ± - ]T ^PxCb+kP . 
The last expression coincides with Eq. (pl6|).B 

6 Conclusions 

We have established the Darboux-covariance of a large class of nonlinear von 
Neumann-type equations. The next step is to employ this fact in construction 
of explicit solutions of such equations. Some classes of solutions have already 
been found in [^0[ for nonlinearities 

ip=[H,f(p)] (53) 

with /(p) = p 2 and f(p) = p q — 2p q ~ 1 (q is an arbitrary real number). Both 
finite- and infinite-dimensional cases were treated by this technique in |2^ ] . 

In a forthcoming paper we will describe other classes of solutions of the 
integrable equations we have introduced. It also seems that Lax pairs that allow 
us to reduce the compatibility conditions to nonlinear equations in a closed form 
can be still generalized. We hope in the future work to develop a description 
of integrable equations of the von Neumann type by taking into consideration 
the Mikhailov method of automorphisms f25j . This type of generalization is 
particularly important if reductions characteristic of Nahm-type equations are 
involved. 
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